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An ab initio analytical model for Fano resonances in two-channel systems is presented. We find 
that the lineshape parameter q factors into a background contribution gb that depends only on the 
uncoupled channels and an interaction contribution gi that is affected by the coupling between the 
channels, revealing how the overall lineshape parameter q may be controlled. In particular, we show 
how conical intersections of the background phase shifts have an important role in the interplay 
between qh and qi . Finally, control of Fano transmission profiles through gb and gi is demonstrated 
for quantum billiards. 
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Fano resonances are ubiquitous whenever discrete lev- 
els interact with a continuum. They have been observed 
in many areas of physics: atomic and molecular systems 
[3, 0] , quantum dots [H , carbon nanotubes [1] , and mi- 
crowave billiards When probed, Fano resonances 
yield their signature scattering cross-section profile [l[ 



which has been used as a fitting function for the results of 
many experiments as well as numerical calculations. Here 
e is the dimensionless reduced energy (which is measured 
in units of the resonance's half-width and shifted so that 
it vanishes at the position of the resonance), while q is 
the Fano lineshape asymmetry parameter. The latter is 
proportional to the ratio between the resonant and non- 
resonant transition amplitudes 0]. Unfortunately, this 
definition does not reveal any means of profile control. In 
this work we suggest how to manipulate the asymmetry 
parameter q through a novel general analytical model for 
two-channels systems that can be specifically applied to 
two dimensional open quantum billiards (QBs). 

Our analysis reproduces the famous result of Eq. Q] 
However, the Fano lineshape parameter q now factors as: 

q = qb-q\- (2) 

Here % is responsible for the contribution of the uncou- 
pled lower (background) channel, whereas qi reflects the 
effect of the interaction between the two channels on the 
resulting profile. Our results show the crucial role played 
by the background transmission unity peaks, which cor- 
respond to conical intersections of the background phase 
shifts, see Fig.[IJa). The factor % corresponds to the ra- 
dial distance from the vertex, and g; behaves like the 
cosine of the azimuthal angle. Finally, the factoriza- 
tion opens the door to resonance profile control which 
we demonstrate for a quantum billiard model. This elu- 




FIG. 1: (a) Parameter space for controlling q near a conical 
intersection of the background phase shifts. The upper cone 
corresponds to the state \E S ) and its phase shift 5s (mea- 
sured relative to \0f). To leading order, this state is an 
azimuth-dependent linear combination of the vertex states 
\Eq^ A ). Changing the coupling V n corresponds to rotat- 
ing the qi "compass wheel", (b) The general form of the 
transmission Fano lineshape. The transmission zero is ob- 
tained at e = — gb9i, while the maximal transmission value 
(1 + gb 2 9i 2 )/(l + gb 2 ) is attained at 1/gbgi. The transmission 
approaches the background value T 1 as |e| — » oo. 

cidates Fano g-reversals in a manner which is more com- 
plete and general than previous attempts Q. 

Two-channel Model. — In order to analyze a multi- 
dimensional scattering process that yields a Fano profile, 
it is often useful to represent the system as a set (albeit 
infinite) of coupled channels. We assume a single-particle 
time-reversal invariant (TRI) multichannel system where 
only the lowest channel is open and the scattering energy 
is in the vicinity of an isolated bound state in the second 
channel, with other bound states energetically separate 
so that we may disregard them. The two-channel Hamil- 
tonian is partitioned as H = H° + V 1 + V n , so that the 
"distorted wave" formalism Q, which treats scattering 
as a conceptual two-stage process, may be applied. The 
first scattering stage, with S-matrix S J , describes scat- 
tering due to V 1 relative to the H° states in the first 
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uncoupled channel. The potential V 1 also accounts for matrix is 
the bound state \ip) at energy E v in the second channel. 
The second scattering stage's S u describes the scattering 
of the H° + V 1 states under the influence of inter-channel 
interaction V 11 , which turns the bound state into a res- 
onance. The scattering matrix for the entire two-stage 
process is S = S I S U @, [!]■ In the present case it is as- ^ 
sumed that S l is known, and we shall show how S* 11 is 
approximated using the generic Fano model 

The first stage is a simple TRI ID problem. Its S- 
matrix may be written as 



,-ti, 



re 
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This matrix has eigenvalues that satisfy e 2lSs — t + r 
and e 2lSA = t — r, hence the notations 5$/a refer to the 
symmetric/antisymmetric combinations of t and r. The 
ambiguity in the above decomposition may be removed 
by one of the discrete gauge choices —tt/2 < 5 r < ir/2 or 
5 s > 5a- We also define the Friedel angle 0p = 5 s + 5a 
the difference AO = 5s — 5a, and the states 



\E») = 



\E A ) = 



V2 

,-iS A 

iV2 



(4a) 
(4b) 



where \E a ) are the first stage scattering states at en- 
ergy E with incoming flux in the direction a = ±. These 
combinations are special in that they have real wavefunc- 
tions with the asymptotic boundary conditions (x\E s ) ~ 
-±= cos(fca; — \5 r ± 5s) and (x\E A ) ~ -j= sin(fcx — 

7}5 r ± 5a) as x — ► ±oo, where the kinetic energy is \k 2 . 
These states are a generalization of the treatment for 
a Hamiltonian with parity symmetry [Io| . The asymp- 
totics also hold for unequal leads by using their respec- 
tive wavenumbers. Time-reversal invariance implies that 
these states' second stage coupling coefficients to the 
bound state v S J, A = {ip\V u \E s / A ) may be chosen to be 
real. Finally, the first-stage transmission coefficient is 



cos - 



(6s -5 A ). 



According to the Fano model [l| , the first stage contin- 
uum is decomposed into states that couple to the bound 
state and those that don't. The coupled states acquire 
the Fano phase shift e 2lSE due to the interaction with the 
bound state, where cot (5e) = —2(E — E v — A E )/T E , 
with the usual definitions T E = 2ir\v E \ 2 and = 
^r^ 2 / dE 1 B ®' E , ■ Here v E is the radial coordinate of the 
point (v E ,v E ) £ R 2 , and we denote its angle as a. 



+iS rn + r I ( 5 ) 



S = S 1 S 11 = e i ( s s+ SA + 5 ®) 1 T 

Ue T '" r /) T r 

t = cos 2 a cos(5s — 5a + 5 E ) + sin 2 a cos(5s — 5a — 6 E ) 
± = cos 2 a sin(<55 — 5 a + 5 E ) + sin 2 a sin(<5s — 5 a — 5 E ) 
±2i sin a cos a sin 5 E . 



If v E varies slowly enough as a function of E then it is 
reasonable to approximate T E and A E using their values 
at the bound state's energy E v These are denoted as 
r and A. In addition, we assume a constant scattering 
background, i.e. S l and its associated parameters are 
(3) treated as constant, with values also taken at £L. This 
is a fair approximation when T is much smaller than the 
scale for which S l changes appreciably. 
-|- f Next, let's look at the Fano lineshapes. The ratio of 

the the full two-stage transmission probability to that of the 
first stage assumes the almost familiar form 



T 

"T 7 



(6) 



with the usual reduced energy e — — cot 5 E and the novel 
factorization of the lincshape parameter in Eq. [2] where 



q h = tan (5 S - 5 A ) = ±y / R I /T I 
q\ = cos(2a). 

In these terms, the actual two-stage transmission coeffi- 
cients is T = (e + <? b <7i)7[(l + e 2 )(l + 9b 2 )]- 

The general form of the transmission Fano lineshape 
is described in Fig. [D(b). Note that a symmetric line- 
shape is obtained not only when the background is zero 
(<7b - * ±oo resulting in a Breit-Wigner Lorentzian) or 
unity (gb = giving a symmetric Lorentzian dip), but 
also when q\ = 0. The latter occurs when the bound 
state is coupled with equal strength to \E S ) and \E A ). 

Conical Intersections. — The lineshape asymmetry pa- 
rameter q factors into a background contribution qb that 
solely depends on V , and a coupling interaction contri- 
bution (/i, which may be controlled using V . However, 
neither q^, nor q^ are gauge invariant when taken alone, 
only their product is. This subtlety will be treated in de- 
tail in what follows. Understanding how V n controls q\ 
is a simple affair — the decoupled first-stage surfaces V 1 
are not affected, so that qi, and the states \E S / A ) remain 
constant regardless of the specific gauge used to define 
them. On the other hand, modifying V 1 may affect both 
gb and gi, an effect that is predominant in the vicinity of 
background unity transmission peaks. 

Suppose that at a scattering energy E = Eq the back- 
ground transmission reaches a unity peak, implying that 
5% = 6 a.' Here we accent values at the transmission peak 



Finally, using these definitions, the total two-stage S- by a so that S\ = e ldp \. Assume throughout that 5s 
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approaches 6 a with a leading order linear in E—Eq. Since 
the background S-matrix is degenerate at this energy, 
there are no preferred eigenvectors. However, we may 
choose to use the states in Eqs. 0] that evolve smoothly 
from those obtained for values E < Eq, which we denote 
\Eq). These limiting states are associated with an an- 
gle 5®, which is just the limit of S r as E is approached 
from below. Note that by a suitable choice of gauge (i.e. 
choosing which point to label x = 0) it is possible to set 
S° r = 0. 

Now we modify V 1 near the background transmission 
unity at Eq, leading to three types of infinitesimal defor- 
mations of the background S-matrix 

dS l = iS\ [ld£ + £° (a x dr) + <r y d()] , 

where £° = exp ( — ia z 6^j. Note the absence of a a z 
generator which breaks TRI. The £ deformation affects 
only 9p, keeping the degeneracy intact, and it doesn't 
change the eigenvectors, nor does it change the position 
of the transmission peak that still reaches unity. The rj 
deformation acts through the a x generator. This removes 
the degeneracy, but to leading order \E^ A ) are still its 
eigenstates. It also shifts the position of the unity peak in 
the transmission spectrum away from Eq. Finally, the £ 
deformations change the height of the peak whilst keep- 
ing its position. More importantly, ^-deformations are 
orthogonal to the //-deformations in the sense that the 
degeneracy is removed so that to leading order the eigen- 
vectors are proportional to \Eq) ± \E A ). 

The effect of these deformations on the first stage S- 
matrix may be described as a conical intersection of the 
phase shifts 8$/a m the T^-plane, schematically depicted 
in Fig. 02a). Here the actual eigenstates \E S I A ) are de- 
fined according to the A9 > gauge, i.e. \E S ) is on the 
upper cone. The leading order behavior of this state de- 
pends on the position in the parameter space. On one 
side of the rj axis \E S ) ~ \Eq), while on the other it 
crosses over to \E S ) ~ \E A ). Similarly, along the £ axis 
it crosses over from \E§) + \E A ) to \E§) ~ \E A ). 

One might naturally ask what topological phase is as- 
sociated with these conical intersections. The answer is 
that 8 r shifts by 2w for each cycle that encircles a ver- 
tex, and this serves as the basis for adiabatic quantum 
swimming/pumping 

The interplay between and qi can now be understood 
in terms of such conical intersections. The A9 > gauge 
choice ensures that </b > on the conical surface (equal- 
ing zero at the vertex), so that up to appropriate param- 
eterization dependent scale factors, q\, corresponds to the 
radial distance from the vertex in parameter space. In 
many cases it is also reasonable to assume that the bound 
state's coupling to the \E^ A ) states remains (relatively) 
constant near the intersection. However, due to the be- 
havior of \E S / A ) that changes according to the position 
in parameter space, the actual couplings to the bound 



state also change, so that q\_ depends on the azimuthal 
angle's cosine. The initial coupling to the \Eq) states 
may be modified through V 11 , and in effect this rotates 
the "compass wheel" in Fig. HJa), i.e. q t = 1 points in a 
different direction in the r^-planc. 

Quantum Billiard Example. — Recently, Fano reso- 
nances have been studied extensively in connection with 
QBs, see e.g. d E3. Usually, the QB is composed of 
an access lead from which the electrons impinge upon 
the main billiard, typically the part associated with a 
quantum dot, and an exit lead. The electron scatter- 
ing problem in QBs can be cast into a one-dimensional 
coupled channel problem 0, [3] , where the channels are 
taken to be the energies of the modes in the transverse 
direction, and the inter-channel interactions are due to 
the non- adiabatic couplings [|, 15|. 

Control of q will be demonstrated on the QB depicted 
in Fig. (Ha). It consists of a rectangular cavity con- 
nected to leads that have two potential barriers of con- 
stant heights Vq, with widths w L / R and distances from the 
cavity l L l R that may be varied to modify V 1 . The lat- 
ter 's first two adiabatic potential surfaces (channels) are 
depicted in Fig.[2jb). The leads may be offset by A L / R in 
order to control the non-adiabatic couplings V 11 result- 
ing from the abrupt change in the local transverse basis 
at the cavity edges. Although the formalism detailed in 
Ref. d is used to calculate the actual non-adiabatic cou- 
plings, symmetry-based selection rules suffice for a qual- 
itative understanding. 
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FIG. 2: (a) Geometry of the quantum billiard with leads used 
to demonstrate control of Fano (/-reversal. Length is mea- 
sured in arbitrary units, while the energy scale is in units of 
[length -2 ] (b) The first two adiabatic surfaces Vl^ix) and the 
second surface bound state \tp). 

The significance of a conical intersection is demon- 
strated near the unity transmission peak at Eq = 27.35 
for a parity symmetric V 1 obtained by setting w L ^ R = 0.1 
and l L l R = 0.295. This is convenient since the symmetry 
ensures that the states \Eq) have respective even and 
odd parity in the S r = gauge. The background V 1 is 
modified to control the lineshape parameter q near the 
conical intersection. Taking A L//fl = 0.01 ensures that V u 
couples \ip) exclusively to \Eq) by symmetry. We choose 
rj = £ L — £ R , which preserves parity, and ^-deformations 
that break this symmetry. 

Figure [3](a) depicts the Fano profiles resulting from 
the numerically exact multichannel calculation [la ], our 
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analytical model, and the background contributions for 
several sets of parameters, see table. In all cases note 
the excellent fit between the actual multichannel calcu- 
lations and the lineshapes predicted by our model. This 
series of calculations corresponds to two g-reversal paths 
in parameter space depicted in Fig. [3Jb) . The first path 
along the r] axis goes from point 1, with q < 0, to 3 
(q > 0) by passing exactly through the background unity 
transmission peak at 2, where q\, = 0. Note that for 
point 1 the transmission peak's energy has moved to the 
right of the resonance, while the opposite happens at 
point 3. The second path goes from 1 to 3 through 2' 
which breaks parity symmetry. This point has a sub- 
unity background transmission peak, corresponding to a 
pure ^-deformation. As expected from the conical inter- 
section, at this point in parameter space the states \E S ^ A ) 
are similar to \Eq ) ± \E^). This implies equal coupling 
to the bound state, so that q\ = 0, giving a symmetric 
dip with the novelty of a less than unity background. 
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FIG. 3: (a) Multichannel and analytical lineshapes obtained 
for various V 1 parameters (see table) . (b) These correspond to 
2 paths in parameter space where q reversal occurs. The over- 
all lineshape parameter vanishes at the intermediate points: 
Qb = at point 2 while qi = at 2'. (c) Variation of the non- 
adiabatic couplings V u , gives g-reversal due to the change in 
the coupling angle a, and hence q\. The table summarizes 
analytical values calculated from the model. 

Control of q may also be affected through V 11 , which 
directly controls q ll whilst qb > remains constant. 
In the language of conical intersections the V 1 param- 
eters remain constant, so that q- reversal is achieved 
through a rotation of the "compass wheel" . This type 
of control works irrespective of the position in parame- 



ter space relative to conical intersections. For simplic- 
ity, again we choose a parity symmetric V 1 , this time 
with i L l R = 0.4125 and w L / R = 0.1 so that q h = 0.544 
throughout. Now we vary A L = —0.03 i— > 0.03, while 
holding A R = 0.03. The results of the numerically exact 
multichannel calculations are depicted in Fig. 03c). The 
results are easily explained in terms of symmetry if we 
use the 8 r = gauge. Since the second channel bound 
state has even parity, it is clear that for A L — —0.03 it 
couples exclusively to \E A ), while for A L = +0.03 it is 
coupled only to \E S ). Equal coupling, with = 0, is 
achieved for A L — 0.00, yielding a sub-unity symmetric 
dip. 

In summary, we have presented an analytical model for 
Fano resonances in coupled two-channel systems, where 
we found that the Fano lineshape parameter q factors into 
background and interaction contributions. The model 
give accurate predictions for the actual transmission line- 
shapes. Moreover, it also provides insight to the relation 
between conical intersections of the background phase 
shifts, the coupling interaction, and the overall lineshape 
parameter q. This allows full control of q, which was 
demonstrated for a quantum billiard example. 
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